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D , We make a further study of CP asymmetries in the pure penguin-induced decay 

modes B' -^ K^*'^^ + K^*'^' and B' -^ K^*^^ + (vr" , p" , a^" ) by using the two-loop 
i^i ■ renormahzation-group-improved effective Hamiltonian and factorization approxima- 

5-H ' tion. Different from the previous results obtained without QCD corrections, the CP 

asymmetries in each of the two groups are enhanced and classified with their respec- 
tive factorization coefficients. A sum over the pure penguin modes of each group is 
proposed to obtain an effective branching ratio and a weighted-mean signal of CP 
violation, which should be statistically significant for experimental observation. 
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1 Introduction 

Observing CP violation in the 5-nieson system and confronting it with the predictions of 
the standard model is an important task in particle physics. The basic signal for CP violation 
is the partial rate difference between a B decay mode and its CP-conjugate process. On the 
basis of the Cabibbo-Kobayashi-Maskawa (CKM) picture, either B^ — B^ mixing [1] or the 
absorptive parts of QCD-loop-induced (penguin) amplitudes [2] can give rise to significant 
CP asymmetries in exclusive nonleptonic B decays. Among various hadronic decay modes 
of B mesons, the pure penguin-induced transitions are of great interest to explore direct 
CP violation in the decay amplitude and to test our understanding of the loop effects 
involved in the quark- level process b ^ q (with q = d or s). Recently evidence for the 
electromagnetic penguin transitions, e.g., B^ — > K*^'y and i?^ -^ K*~'y, has been obtained 
by CLEO Collaboration [3] . Accordingly studying CP violation in the pure strong penguin- 
induced B decays becomes a necessary and realistic topic today. 

In the literature [4-6], CP asymmetries in some pure penguin channels have been calcu- 
lated with the Q CD- uncorrected effective Hamiltonian H^^^ (see Eq.(2)) and factorization 
approximation. For the processes occurring through b —^ d (or b ^ s), those results are 
approximately independent of the spin properties of final-state mesons. Hence a sum over 
the relevant modes has been proposed in order to reduce the total number of B events and 
to obtain a statistically significant signal of CP violation. In this paper, we shall use the 
two- loop renormalization-group-improved effective Hamiltonian TYcff (see Eq.(4)) to make a 
further study of CP asymmetries in the following two handfuls of pure penguin transitions: 

(I) B- — > K^K-, K*°K-, K^K*-, K*^K*- ; 

(II) B- — > K'^n-, K^p-, K°ai, K*^tx" , K*^ p' , K*^a^ . 

It is worth emphasizing that the above ten decay modes are flavor self-tagging processes, 
which should be favored for experimental reconstructions. In addition, there are no QCD- 
loop-induced hairpin diagrams contributing to these transition^ and the electroweak penguin 
effects on them are negligible [8]. Thus application of the factorization approximation to 
these decays should lead to less uncertainty. Different from the naive results obtained with 
H^g^, the magnitude of CP asymmetries in the transitions (I) and (II) are enhanced by the 
next-to-leading order QCD corrections from Ties- In each group the asymmetries become 
non-degenerate and are classifled with their respective factorization coefficients, but they 

^In our sense, some of the QCD-loop-induced transitions such as B^ -^ K^cf) are not the pure penguin 
decay modes. The reason is that they can also occur through the hairpin diagrams. For a detailed discussion, 
see Ref. [7]. 



still keep the same sign. Hence we propose a sum over the modes of each group to obtain 
an effective branching ratio and a weighted-mean signal of CP violation, which should be 
statistically significant for observation. The feasibility of measuring such weighted-mean CP 
asymmetries the uncertainties in evaluating them are briefly discussed. 



2 Effective Hamiltonians and factorization approximation 

To calculate decay amplitudes of the aforelisted pure penguin channels, we use the low- 
energy effective Hamiltonians for Ai? = ±1 transitions and the well-known factorization 
approximation [9]. Neglecting influence of the renormalization group, one used to apply the 
one-loop penguin Hamiltonian [4-6] 
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to phenomenological discussions. In H^^'\ the QCD renormalization effect is assumed to be 
approximately included by the effective coupling constant cts at the physical scale /i = mj,; 
Vi = VibV*^ (with i = u,c,t and q = d,s) are the CKM factors corresponding to 6 ^ g; FiihP') 
are the loop integral functions of the time-like penguin diagrams versus the virtual gluon 
momentum A;^ [2,6]; Nc is the number of colors; and (53,...,6 represent the penguin operators. 
Together with the current-current operators Qi 25 Q^,...fi form an operator basis as follows: 

Qi = {QaUp)y_A{uph^)y_A , Q2 = {qu)v-A{ub)v-A , 

Qi = {QaC^)v-A{C(}ba)v-A > Q2 = {Qc)v-A{cb)v-A , 

Qs = {qb)y_AY.i^'l')v-A, Qi = iQabf3)v-AY.i^'(3l'a)v-A , (3) 
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Compared with H^g'^, the two-loop QCD-corrected Hamiltonian Ti^s is of the form [10] 

Heff (Afi = -1) = ^ Vu fe C.Qr) + Vc fe C^Q?) - Vt fe ^^Q^] " (4) 

Here the Wilson coefficient functions Ci^,,,^q, obtained by applying the renormalization-group- 
improved perturbation theory, have included the next-to-leading order QCD corrections at 
the scale fi = mi,. Note that q depend on the renormalization scheme used for the four- 
quark operators Qi. This scheme dependence, however, can be cancelled by certain one-loop 
matrix elements of the decay mode in question. For a detailed description of the approach to 



obtain the renormalization-scheme independent transition amplitudes in two-body B decays, 
we refer the reader to the article of Fleischer [11]. 

With the help of Hcs, the amplitude of a two-body pure penguin- induced decay mode 
B^ -^ XY may be expressed as linear combinations of < XY\Qi\B^ > with weighting 
factors Cj and fj. Following Ref. [11], one can remove the renormalization-scheme dependence 
of Cj in the physical transition amplitudes by taking one-loop penguin matrix elements of the 
operators (^2''^ into account (see Fig. 1). The renormalization-scheme independent amplitude 
of B^ -^ XY is finally obtained as [11,12] 
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where Cj (with i = 2, ..., 6) are the scheme independent Wilson coefficients including the next- 
to- leading order QCD corrections [10], and Mp and Mp are the hadronic matrix elements of 
penguin operators defined by 
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In contrast, the decay amplitude of B^ -^ XY calculated with H"^^^ is of the form 
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The hadronic matrix elements Mp and Mp are calculated using the naive factorization 
approximation [9]. Here we only consider contributions from the spectator- type diagrams 
as illustrated in Fig. 1, while the OZI- forbidden transitions and (or) formfactor-suppressed 
annihilation topologies are neglected. Consequently Mp and Mp are both dominated by a 
single term 
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and their factorized expressions are given as 
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In Eq. (9), the factorization coefficient C,p is obtained by transforming the {V — A)(y + A) (or 
{S + P){S — P)) currents of Qs.e into the (V — A){y — A) ones. For each specific decay mode. 



(p depends upon the current masses of quarks as well as the angular momenta and parities 
of X and Y mesons. We list the expressions of (p in Table 1 for 0^0^, 0^1''^, 1^0^, and 1^1''= 
final states, respectively. The hadronic matrix elements MSC, can be Lorentz-invariantly 
decomposed with the relevant decay constants and formfactors, as given in Refs. [6,12]. For 
our purpose, we shall sum over the polarizations of vector or axial vector mesons in the final 
state and estimate |M^yP by means of the Bauer-Stech-Wirbel (BSW) model [13]. 



3 CP asyminetries 

Now let us define the CP-violating partial decay rate asymmetry for B~ -^ f and its 
CP-conjugate counterpart P^ -^ f: 

, . f. _ r(a+ ^ /) - r(B- -* /) 

In the CKM mechanism with three families of quarks, every of the aforementioned pure 
penguin transitions contains two amplitude components that have different CP-violating 
weak phases (f„ and Vc) and different CP-conserving strong phases (ImPu(A;^) and ImPc(^^))- 
The CP asymmetry is approximately independent of hadronic matrix elements, since M^,^, 
can be cancelled in Acp{f)- Thus Acp{f) is simply expressed as [6] 
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where P/ is a ratio of strong-interaction parts of the two amplitude components correspond- 
ing to Vc and t>„. When calculating Acp{f) with P^g°, we obtain [4-6] 

' ~ F^{k^) - P,(P) • ^''> 

This implies that CP asymmetries in the transitions (I) or (II) will have the same magnitude 
and the same sign. While applying Hes to the analysis, P/ becomes 

, _ Fcje) - sjcp) 
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Clearly Pj includes the next-to-leading order QCD corrections and depends upon the fac- 
torization parameter (p. Since the values of (p for B^ -^ K^K^ ,K^K*^, and K*^K^ (or 



K*^K*~) are different from one another, their CP asymmetries will become non-degenerate 
in contrast with the result obtained from Eq. (12). So is the situation for the pure penguin 
channels B^ -^ ^(*)o -|- (tt^, p^, ai). This feature, arising from the QCD corrections, requirs 
us to reconsider the idea of suming over the pure penguin transitions (I) or (II) to gain a 
statistically significant CP-violating signal. 

A comparison between R^ and Rf shows that the QCD improvement only modifies 
R-e-Fu,c(^^)- Hence the asymmetries in each group of decay modes remain the same sign. 
This ensures that there should be little dilution effect on the CP asymmetries if one sums 
over the available final states / and / in the transitions (I) or (II). In this case, one can 
obtain an effective branching ratio like 
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where Tq is the lifetime of B^ meson. The corresponding CP asymmetry Acp is a weighted- 
mean value of Acp{f)'- 
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Although the value of Acp is impossible to deviate too much from the ones of Acp{f), it 
is expected that the effective branching ratio B^s can become several times larger than a 
single B{f). Thus the total number of B^ events needed in measuring Acp will be remark- 
ably reduced and could be accessible in the forthcoming B-meson factories. Note that the 
weighted-mean asymmetry depends upon hadronic matrix elements and is therefore difficult 
to be evaluated in a reliable way. In the long run, however, a great improvement of the 
present calculations will be possible to yield trustworthy results. 



IV Numerical results and discussion 



For illustration, we make a numerical estimate of CP asymmetries and branching ratios 
for the pure penguin-induced decay modes (I) and (II). The CKM parameters are taken as 
A = 0.22, A = 0.90, p = —0.50, and r] = 0.30, which are consistent with the present data on 
B^ — B^ mixing and e^ [14]. We choose the values of the current quark masses as m„ = 5 
MeV, md = 10 MeV, m, = 175 MeV, m^ = 1.35 GeV, nib = 4.8 GeV, and rrit = 150 GeV. 

(5) 

Fixing A^^ = 200 MeV, the effective QCD coupling constant and next-to-leading order 
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Wilson coefficients are obtained [10] as as{mb) = 0.18, ci = —0.273, C2 = 1.123, C3 = 0.014 
C4 = —0.032, C5 = 0.009, and cg = —0.039. With the approximate rule of discarding 1/A^, 
corrections in the naive factorization approach [9], here we take 1/Nc = for the relevant 
charmless nonleptonic B transitions. It is not clear what value of k'^ (the virtual gluon 
momentum in penguin graphs) should be taken in exclusive B decays. A simple estimate, 
involving two-body kinematics and one-gluon exchange to accelerate the spectator quark, 
yields ml/4 < k"^ < vn^jl [4]. To evaluate branching ratios we use r^ = 1.29 x 10^^^ s [15] 
and quote values of the relevant decay constants and formfactors from Ref. [13]. 

Our numerical results of CP asymmetries Acpi^]^ and branching ratios -B(/) are given 
in Fig. 2 and Tables 2 and 3, respectively, as functions of k'^ . The following features can be 
observed: 

(1) Compared with the previous result without the next-to-leading order QCD correc- 
tions, now CP asymmetries in the transitions (I) or (II) are enhanced to some extent and 
become non-degenerate for the 0^0^, 0^1^, and 1^0^ (or 1^1^) final states. As expected, 
-^cp(/) keep the same sign even though their values may change significantly with k^ . 

(2) Among the transitions (I), the 0^1^-type decay mode B^ -^ K^K*^ has the largest 
CP asymmetry but the smallest branching ratio. In contrast, the 0^0^-type channel B^ -^ 
K^K~ has a relatively small CP asymmetry in spite of its largest branching ratio. So is the 
situation for the 0^1"^- and 0^0^-type decay modes in the transitions (II). One can find that 
there exists no individual channel which has obvious advantages over the others of the group 
for measuring CP violation. For this reason, studying the weighted-mean CP asymmetries 
should be interesting in statistics. 

(3) For either the transitions (I) or (II), the weighted- mean asymmetry Acp is almost 
of the same magnitude as each Acp{f), but its corresponding effective branching ratio 
Bcs is several times larger than the individual B{f). For example, B^q/ B{K^K^) ~ 2.5, 
B^s/B{K^K*-) ~ 7, B^s/B{K*^K~) ~ 4, and B^^/ B{K*^K*-) ~ 3.5. As a result, mea- 
suring Acp instead of Acp{f) should reduce the needed B^ events several times and is 
therefore favored in practice. At the 3o" level, our estimates indicate that about 5 x 10^ 
and 2 x 10^ B^ pairs are available to probe the weighted-mean CP asymmetries of 7% and 
1%, respectively, in the transitions (I) and (II). In experiments, to detect a handful of pure 
penguin channels with similar final-state particles (e.g., B~ -^ f^(*)o _|_ f^(*)-) might not be 
more difficult than to detect an individual channel among them. 

Except for the uncertainty induced by the CKM parameters, there are some other theo- 
retical uncertainties from the assumptions and approximations made in our calculations. It 



is worthwhile at this point to give a brief summary of them: 

(1) For the pure penguin-induced exclusive B decays in question, CP violation has been 
treated by postulating that the phases of the penguin amplitudes at the meson level is the 
same as those of the penguin loops at the quark level. This leads to the problem that the CP 
asymmetries Acp{f) depend strongly upon the virtual gluon momentum (A;^) of the timelike 
penguin graphs. The large uncertainty arised from A;^ cannot be experimentally limited in 
a meaningful way, since it is not a measurable. To modify this unsatisfactory point, an 
attempt has been made by calculating the k'^ distribution of some exclusive B decays and 
fold it with the momentum dependence of the loop amplitudes [16]. 

(2) Neglecting inelastic rescattering effects, pure penguin decay modes have only one 
isospin component. Hence their CP asymmetries should suffer little from the unknown 
final-state interactions. Under the same condition, a sum over the transitions (I) or (II) is 
in principle available without suffering cancellations. Whether the above simplification is 
reasonable or not can be clarified in the forthcoming experiments of 5-meson physics. 

(3) It is seen that the CP asymmetries Acp{f) are independent of hadronic matrix 
elements in the factorization approximation. However, the branching ratios B{f) should be 
sensitive to the model applied to hadronic matrix elements. Although the weighted-mean 
asymmetries Acp are also affected by hadronic physics, the involved uncertainty should not 
be drastic in general. The reason is that both Acp{f) and Acp are ratios of the partial 
decay rates. As an illustration, our numerical estimates of Acp and Bcs could give one a 
feeling of ballpark numbers to be expected. 

In conclusion, the pure penguin-induced B^ decays are of great interest to uncover direct 
CP violation in the decay amplitude. A weighted-mean signal of CP asymmetries may be 
more easily observed in experiments, although to evaluate it with reliability remains difficult 
in theories. In order to test the strong-interaction penguin picture and direct CP violation 
mechanism at 5-meson factories, a further study of the dynamics of nonleptonic exclusive 
B decays is urgent. 
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Figure Captions 

Fig. 1 Pure penguin-induced decay modes 5^ -^ X{qq') + Y{q'u) with (g, q') = (rf, s) 
or (s, d) : (a) the tree-level matrix elements of the penguin operators, and (b) the one-loop 
penguin matrix elements of the current-current operators. 

Fig. 2 CP asymmetries versus k'^ (virtual gluon momentum in penguin graphs) for a 
few pure penguin- induced B^ decay modes: (a) exclusive b —^ d, and (b) exclusive 6 — > s. 
The circled curve is the degenerate result without QCD corrections; the dashed, dotted, 
and dot-dashed curves stand for the QCD-corrected results; and the dagged curve is the 
weighted- mean CP- violating signal. 



Table Captions 

Table 1 Factorization parameters (p (defined in Eq. (9)) for tlie decay modes 5^ -^ 
XY, wliere X and Y may be 0"*^ or 1"*= mesons. 

Table 2 Brancliing ratios B{f) (xlO^^) versus k'^ (virtual gluon momentum in penguin 
graphs) for B^ -^ 7^(*)o -|- _ft'(*)^. B^s (xlO^'^) is an effective branching ratio obtained by 
summing over these four pure penguin modes. 

Table 3 Branching ratios B{f) (xlO~®) versus k"^ (virtual gluon momentum in penguin 
graphs) for B^ -^ ^(*)o -|- (tt", p~, a{). Bcs (x 10^^) is an effective branching ratio obtained 
by summing over these six pure penguin modes. 
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Table 2 



k^ml 


B{K^K-) 


B{K^K*-) 


B{K*^K-) 


B{K*^K*-) 


Bcs 


0.1 


12.5 


0.39 


7.47 


9.20 


29.6 


0.2 


12.5 


0.39 


7.47 


9.20 


29.6 


0.3 


14.2 


0.46 


8.54 


10.5 


33.7 
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15.8 


0.53 


9.56 


11.8 


37.7 
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14.6 


0.49 
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34.8 
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0.46 
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13.1 
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31.1 
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12.5 


0.41 


7.53 


9.27 


29.7 
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12.0 


0.39 


7.22 
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28.5 


1.0 


11.6 


0.38 


6.98 
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27.6 
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0.1 


7.77 


0.22 


0.19 


4.84 


5.99 


5.39 


24.4 


0.2 


8.77 


0.25 


0.23 


5.50 


6.81 


6.13 


27.7 


0.3 


11.3 


0.36 


0.32 


7.20 


8.92 


8.03 


36.1 


0.4 
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0.39 


0.35 


7.67 


9.50 


8.55 


38.5 


0.5 


10.7 


0.34 


0.30 


0.81 


8.43 


7.59 


34.2 


0.6 


9.79 


0.31 


0.28 


6.23 


7.71 
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31.3 
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0.28 


0.25 
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7.17 
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29.1 
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0.27 


0.24 
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27.5 
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0.22 
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6.43 


5.79 


26.1 


1.0 


7.90 


0.24 


0.21 


4.99 


6.17 
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25.1 



